arXiv: 1501.0391 lv2 [nucl-th] 30 Sep 2015 


The structural evolution in transitional nuclei of mass 82 < A < 132 

M. Bhuyan* 

State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 

(Dated: October 1, 2015) 

In this theoretical study, we report an investigation on the behavior of two-neutron separation energy, a differ¬ 
ential variation of the nucleon separation energy, the nuclear charge radii and the single-particle energy levels 
along the isotopic chains of transitional nuclei. We have used the relativistic mean field formalism with NL3 
and NL3* forces for this present analysis. The study refers to the even-even nuclei such as Zr, Mo, Ru and Pd 
for N = 42—86, where a rich collective phenomena such as proton radioactivity, cluster or nucleus radioactivity, 
exotic shapes. Island of Inversion and etc. are observed. We found that there are few non-monotonic aspects 
over the isotopic chain, which are correlated with the structural properties like shell/sub-shell closures, the shape 
transition, clustering and magicity etc. In addition to these, we have shown the internal configuration of these 
nuclei to get a further insight into the reason for these discrepancies. 


PACS numbers: 21.10.Dr, 21.10.Ft, 21.10.Gv, 21.10.Tg 

I. INTRODUCTION 

Nowadays, one of the most sensitive and crucial region in 
the nuclear chart for investigation is laying in between Z= 
35—64 and A= 82—132. This region reveals a large num¬ 
ber of interesting discoveries of new phenomena, such as 
proton radioactivity 0]^, cluster radioactivity HUH, exotic 
shapes EH, Island of inversion iiiioi, abnormal varia¬ 
tion of major shell closures (i.e. extra stability near drip¬ 
line) lEHia and giant halo near neutron drip-line region m 
etc. These crucial features may be due to the rapid grow¬ 
ing possibility of the neutron-proton ratio {N/Z) in a nu¬ 
cleus. From last few decades, it is possible to study these 
exotic nuclei by using the radioactive isotope beams (RIB) fa¬ 
cilities. This reveals the new concept entitled as aforemen¬ 
tioned magic number. In other word, the confirmation of 
magic number near /3—stability line are not mandatory uni¬ 
versal iUMni- Further, the structural properties of nuclei 
far away from the /3—stability line are also active areas of 
research in both theories and experiments ifTJl fT4l flSl . In 
particular, the neutron-rich Zr—, Mo—, Ru— and Pd— with 
mass numbers A =100—130 are of special interest for various 
reasons. For example, they lie far away from the /3—stable 
region of Nuclear Landscape, result in a well established de¬ 
formation, but close enough in the magnitude of microscopic 
excitations to compete with the collectivity of double shell 
closure nuclei IHIIllIlQl- Moreover, these nuclei are also 
holding an active participation in the nucleosynthesis of heavy 
nuclei in astrophysical r-process. The mass and decay prop¬ 
erties are quite an essential ingredient to building up the path, 
the isotopic abundances and the time period of these process 

ED. 

In addition to that the nuclear structure of these nuclei is 
characterized by a strong competition between various shapes, 
which gives rise to the shape instabilities that lead to coex¬ 
istence nuclear shape transitions in the isotopic chains fTfX . 
This could be understood from the potential energy surface 
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at different deformations. Elaborately, the occurrence of two 
(or more) nearly equally deep minima in the potential energy 
surface at different deformations show the signature for nu¬ 
clear shape coexistence. Hence, one can say the nuclear shape 
are not only vary with the nucleon number but also with the 
excitation energy and spin. It is well known that the bind¬ 
ing energy of a nucleus is one of the most precise measured 
observable from the experiments ll23l l24ll . Several nuclear 
observables which are highly relevant for understanding var¬ 
ious features of the nuclear structure can be computed from 
its mass such as the average nuclear field, nucleon-nucleon 
(NN) potential, single particle energy etc. The correlations 
among these fundamental quantities are amended to explain 
the deformed ground states, low-lying isomeric states and few 
derived quantities like moments of inertia and vibrational ex¬ 
citation energy etc Il25l428l . It is acclaimed that the energy 
involved in removal of Fermions from a strongly correlated 
system of identical Fermions must be a good indicator for the 
stability of the system. This magnitude of this energy has 
much higher values for systems with even number of parti¬ 
cles than odd one, if the pairing is a dominant component in 
the binary Fermion—Fermion interaction. 

In this present work, the quantities of interest are the nu¬ 
clear potential energy surface, nuclear shape, nuclear binding 
energy, two neutron separation energies (S 2 n), the differential 
variation of neutron separation energy AS 2 n, the root-mean- 
square charge distribution rch and the single particle energy 
level for the even-even mass transition nuclei. Base on these 
decisive observables, we have focused on the evolution on the 
structural properties of transition nuclei. The paper is orga¬ 
nized as follows: Section II gives a brief description of the 
relativistic mean field formalism including the pairing energy 
correlation. The results of our calculation along with discus¬ 
sions are presented in Section III. Section IV includes a short 
summary along with few concluding remarks. 


II. THE RELATIVISTIC MEAN-FIELD (RMF) METHOD 

From last few decades, the nuclear covariant density func¬ 
tional theories (CDFT) are quite successful in describing the 
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ground and the intrinsic excited state including fission states 
of the exotic heavy and superheavy nuclei over the nuclear 
chart ||29143^ . Basically, there are four different patterns to 
perform the covariant density functional: the point coupling 
nucleon or meson exchange interactions connected with the 
density-dependent or nonlinear couplings. One can also in¬ 
troduced all shape degrees of freedom to CDFTs by breaking 
both the axial and reflection symmetries simultaneously (see 
the Ref. 041 for more details). In the relativistic mean field 
approach, the nucleus is considered as a composite system of 
nucleons (proton and neutron) interacting through exchanges 
of mesons and photons OTIOTWSI . Further, the contributions 
from the meson fields are described either by mean fields or 
by point-like interactions between the nucleons ll43l |4^ . and 
the density dependent coupling constants Il45ll46l or nonlinear 
coupling terms ll47l |48l are introduced to reproduced the cor¬ 
rect saturation properties of infinite nuclear matter. Flere, most 
of the computational efforts are devoted to solving the Dirac 
equation and to calculate various densities. In the present cal¬ 
culation, we have used the microscopic self-consistent rela¬ 
tivistic mean field (RMF) theory as a standard tool to investi¬ 
gate the nuclear structure phenomena. It is worth mentioning 
that the RMF approach is one of the most popular and widely 
used formalism among them. The relativistic Lagrangian den¬ 
sity (after several modification of the original Walecka La¬ 
grangian to take care of various limitations) for a nucleon- 
meson many body system EMISlHSa is given as; 

C = - M}ipi + 

+ lrnlV^^V^ - 

( 1 ) 

From the above Lagrangian, we obtain the field equations for 
the nucleons and mesons. These equations are solved by ex¬ 
panding the upper and lower components of the Dirac spinors 
and the boson fields in an axially deformed harmonic oscil¬ 
lator basis, with an initial deformation Pq. The set of cou¬ 
pled equations are solved numerically by a self-consistent it¬ 
eration method HUHST). Based on the effective interactions 
used in the RMF functional, the center of mass energy can 
be calculated either in harmonic oscillator approximation or 
from the quasi-particle vacuum self-consistently. In case of 
oscillator approximation, the spurious center of mass motion 
is substracted using the Elliott-Skyrme approximation 1581 . 
The analytical form is given as: 

E,,^ = ^A-K (2) 

where A is the mass number. In other hand, one should esti¬ 
mate the center of mass energy using self-consistent method 
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FIG. 1: (Color online) The obtained results for nuclear relative bind¬ 
ing energy difference {AE)b = Eb{Nf = 18) — Eb{Nf = 
10 — 22), charge radius Teh and the quadruple deformation P 2 for 

100,1262^, 102>128Mo, 106.132p^ 

and NL3* force parameter. See text for more details. 
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where, \E) = \E)jimf wave function. The P and A are the 
total linear momentum and the nuclear mass number, respec¬ 
tively. The results obtained from these two methods for Zr 
isotopes are given in Table 1. From the table it is clear that, 
the calculated center of mass energies from both the cases are 
almost overlap with each other. Hence, one can use any one 
of the method for the center of mass energy correction in the 
calculation of this nuclear mass region (see Ref. iOl). The 
total quadrupole deformation parameter P 2 is evaluated from 
the resulting proton and neutron quadrupole moments, as 


Q = Qn+Qp = ^J^{^AR^P2). (4) 

The root mean square (rms) matter radius is defined as 

(rV) = j J p{rj_,zydT, (5) 

where A is the mass number, and p{r±, z) is the axially de¬ 
formed density. The total binding energy and other observ¬ 
ables are also obtained by using the standard relations, given 
in B9|| . As outputs, we obtain different potentials, densities, 
single-particle energy levels, radii, deformations and the bind¬ 
ing energies. For a given nucleus, the maximum binding en¬ 
ergy corresponds to the ground state and other solutions are 
obtained as various excited intrinsic states. 
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TABLE I: The center-of-mass energy obtained from phenomenological l58l and microscopic self-consistent methods I59II60I for Zr isotopes. 


Method 

loozr 

i02zr 

104z, 

i06zr 

lOSzr 

“°Zr 


“"Zr 

“®Zr 

“»Zr 

i20zr 

Ec.m = |A"3 

-6.63 

-6.58 

-6.54 

-6.46 

-6.46 

-6.42 

-6.38 

-6.34 

-6.31 

-6.27 

-6.23 

p _ iF\P^\F) 

-C/c.m — 2M 

-6.83 

-6.79 

-6.75 

-6.68 

-6.61 

-6.55 

-6.49 

-6.45 

-6.41 

-6.44 

-6.47 


A. Pairing Energy 


To deal the nuclear bulk properties of open-shell nuclei, one 
can not neglect the pairing correlation in their ground and 
intrinsic excited state IHl. There are various methods such 
as the BCS approach, the Bogoliubov transformation and the 
particle number conserving methods have been developed to 
treat the pairing effects in the study of nuclear properties in¬ 
cluding fission barriers Il62l - l65l . In principle, the Bogoliubov 
transformation is widely used methods to take pairing corre¬ 
lation into account for the drip-line region Il39l - l42l 1^ l67l . In 
case of nuclei not too far from the /3-stability line, one can use 
the constant gap BCS pairing approach reasonably good to 
take care of pairing ll68]l . Further, the BCS approach may fail 
for neutron-rich light nuclei. In the present analysis, we have 
considered the intermediate mass neutron-rich nuclei, hence 
the RMF results with BCS treatment should be reliable. In 
other word, to avoid the difficulty in the calculation, one can 
employ the constant gap BCS approach to deal the present 
mass region EIllhaiTOl. Now the expression for pairing en¬ 
ergy is given by. 


^pair — G 


E 

.i>0 


UiVi 


( 6 ) 


where vf + = 1 is known as the occupation probability 

and G is the pairing force constant ll7TI473l . The variational 
procedure with respect to the occupation numbers v^, gives 
the BCS equation 2eiUiVi — A(m| — vf) =0 and the pairing 
gap A is defined by 


G is not calculated explicitly from the RMF equations. Using 
the above gap parameter, we calculate the occupation proba¬ 
bility and the chemical potentials A„ and Xp from the particle 
numbers using the above equations. Now, we can rewrite the 
pairing energy as, 

A2 ^ 

Epair — — A y ^ UiVi- (9) 

z>0 

Since it depends on the occupation probabilities vf and uf, 
the pairing energy should change with particle number for a 
constant pairing gap. It is well known that the pairing en¬ 
ergy Epair diverges if it is extended to an infinite configu¬ 
ration space for a constant pairing gap A and force constant 
G. Also, for the states spherical or deformed, with large mo¬ 
menta near the Fermi surface, A decreases in all the realistic 
calculations with finite range forces. However, for the sake 
of simplicity of the calculation, we have assumed, the pair¬ 
ing gap for all states is equal near the Fermi surface. In the 
present calculations we have used a pairing window, and all 
the equations extended up to the level — A < 
where a factor of 2 has been included in order to reproduce the 
pairing correlation energy for neutrons in ^^®Sn using Gogny 
force CD El. This kind of approach has already been used by 
many other authors in RMF and Skyrme-Hartree-Fock (SHF) 
models EU EH El El- 


III. DETAILS OE CALCULATION AND RESULTS 
DISCUSSION 


A = gY,W ( 7 ) 

i>0 


This is the famous BCS equation for pairing energy. The den¬ 
sities are contained within the occupation number. 


Ui = V, = 


1 - 


Ci - A 


y/(e.-A)2 + A2 


( 8 ) 


The standard expression for the pairing gaps of proton 
and neutron are Ap = and A„ = 

respectively 16^ . Here the constants and 
their values are as follows: R = 5.72, s = 0.118, f = 8.12, i?s = 
l,and/ = {N—Z)/{N + Z). (Note that the gaps obtained by 
these expressions are valid for nuclei both on or away from the 
stability line for this mass region). The pairing force constant 


In the present work, we have used the successful NL3 ll74ll 
and the recently proposed NL3* ITTSII force parameters, which 
are excellent in the description of ground and excited states 
with many collective aspect for spherical and deformed nu¬ 
clei. In the mean time, there are several other mean-field inter¬ 
actions have been developed. In particular, the density depen¬ 
dent meson-exchange DD-MEl ll46l and DD-ME2 fT^ in¬ 
teractions, which are adjusted to improve the isovector chan¬ 
nel. Eurther, the density dependent points coupling interaction 
01 El has been developed to describe the deformed heavy 
and superheavy nuclei. Even these interactions have been de¬ 
veloped to provide a very successful description of various 
special features. At present, the AL3 ll74ll and NL3* ITTSII 
forces are also accepted in compete with these parameters 
to reproduce the properties of the stable and nuclei far from 
the /3-stability line. In RMET, the mean-field equations are 
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TABLE II: The binding energy Eb, root-mean-square charge radius Vch and the quadrupole deformation parameter P 2 for the ground states 
and few selective first intrinsic ecxited state of ® 2 -i 262 j. nuclei compare with the experimental data 1821 . wherever available. 

See the text for more details. 


N 

RMF (NL3) 

RMF (NL3*) 

Experiment 

RMF (NL3) 

RMF (NL3*) 

Experiment 


BE 

f ch 

02 

BE Teh 

P 2 

BE rch 

02 

BE 

^ ch 

02 

BE 

rch P 2 

BE rch 

P 2 

82-126^j. 

86-130jjy 

42 

691.3 

4.29 

-0.197 

690.9 4.29 

-0.192 


0.367 

698.7 

4.39 

-0.204 

698.1 

4.39 -0.199 




691.0 

4.28 

0.488 

690.6 4.28 

0.488 










44 

715.7 

4.29 

-0.188 

715.8 4.27 

-0.001 


0.251 

726.8 

4.39 

0.053 

726.8 

4.38 0.058 




715.5 

4.28 

0.469 

715.5 4.28 

0.467 










46 

739.2 

4.28 

0.001 

739.2 4.27 

0.001 

740.6 

0.151 

755.7 

4.39 

0.096 

755.5 

4.38 0.089 



48 

762.6 

4.28 

0.001 

761.9 4.28 

0.001 

762.6 4.2812 

0.185 

782.3 

4.37 

0.006 

782.0 

4.37 0.004 



50 

783.9 

4.28 

0.001 

783.1 4.28 

0.000 

783.9 4.2696 

0.089 

808.0 

4.37 

0.001 

807.4 

4.37 0.000 

806.9 


52 

797.8 

4.29 

0.001 

797.2 4.28 

0.000 

799.8 4.3057 

0.1027 

825.6 

4.39 

0.003 

825.1 

4.38 0.002 

826.5 4.3927 

0.1579 

54 

810.5 

4.34 

0.169 

808.9 4.31 

0.002 

814.7 4.3312 

0.090 

843.2 

4.42 

0.159 

842.5 

4.42 0.155 

844.8 4.4232 

0.1947 

56 

824.5 

4.38 

0.243 

822.9 4.38 

0.233 

828.9 4.3498 

0.080 

860.1 

4.45 

0.205 

859.2 

4.44 0.199 

861.9 4.4536 

0.2148 

58 

837.0 

4.42 

0.318 

834.9 4.40 

0.274 

840.9 4.4185 


875.3 

4.47 

0.225 

874.3 

4.47 0.216 

877.9 4.4818 

0.2404 

60 

849.8 

4.48 

0.432 

847.6 4.49 

0.453 

852.2 4.5220 

0.355 

889.3 

4.49 

0.234 

888.2 

4.48 0.215 

893.0 4.5104 

0.2707 

62 

860.6 

4.50 

0.428 

858.2 4.50 

0.428 

863.7 4.5690 

0.427 

903.9 

4.57 

0.385 

901.7 

4.52 0.295 

907.5 

0.257 

64 

870.6 

4.52 

0.427 

868.0 4.52 

0.424 

873.8 

0.38 

917.6 

4.54 

-0.236 

916.3 

4.53 -0.232 

920.9 

0.292 









916.9 

4.58 

0.373 

915.8 

4.57 0.371 



66 

880.4 

4.54 

0.419 

877.6 4.54 

0.418 

883.9 


931.0 

4.55 

-0.236 

929.5 

4.55 -0.234 

933.3 

0.295 









929.1 

4.59 

0.357 

928.0 

4.59 0.357 



68 

889.8 

4.56 

0.416 

886.8 4.56 

0.419 

892.6 


943.4 

4.57 

-0.239 

941.5 

4.57 -0.238 

945.0 

0.306 









941.5 

4.60 

0.349 

940.1 

4.60 0.348 



70 

897.2 

4.59 

0.445 

893.9 4.59 

0.461 

900.4 


954.1 

4.59 

-0.233 

951.8 

4.58 -0.229 



72 

904.1 

4.62 

0.478 

900.4 4.62 

0.479 



963.9 

4.59 

-0.203 

961.5 

4.59 -0.196 



74 

911.8 

4.52 

-0.170 

908.8 4.52 

-0.166 



974.2 

4.61 

-0.178 

971.5 

4.60 -0.176 



76 

917.7 

4.52 

-0.109 

914.5 4.52 

-0.095 



983.1 

4.62 

-0.169 

979.9 

4.61 -0.167 



78 

923.4 

4.52 

0.065 

920.0 4.52 

0.043 



991.6 

4.61 

0.114 

988.7 

4.60 0.111 



80 

929.6 

4.54 

0.003 

925.7 4.53 

0.002 



999.5 

4.61 

0.074 

996.2 

4.61 0.074 



82 

935.6 

4.55 

0.001 

932.0 4.56 

0.001 



1007.3 

4.62 

0.001 

1003.4 

4.61 0.001 



84 

936.7 

4.56 

0.001 

932.8 4.57 

0.064 



1010.1 

4.64 

0.067 

1006.3 

4.64 0.076 



86 

937.3 

4.57 

0.069 

933.6 4.58 

0.079 



1013.4 4.66 

0.139 

1009.4 4.66 0.138 




solved self-consistently by taking different inputs of the ini¬ 
tial deformation called Po igEiiiTiiiTiiiTa. To verify the 
convergence of the ground state solutions for this mass re¬ 
gion, we pursued the calculation for Ng =20 and varying Np 
from 10 to 22. The difference between the binding energy ob¬ 
tained from Np=12i to Nf=IQ-22 is entitled as relative bind¬ 
ing energy difference and denoted as {AE)b- The estimated 
relative binding energy difference {AE)b = Eb{Np = 
18) — Eb{Np = 10 — 22), the charge radius Vch and the 
quadruple deformation P 2 for ^ 00 , 1262 ^^.^ i 02 .i 28 ]yjQ I04,i30]^jj 
and i06,i32pj isotopes from NL3 and NL3* force are shown 
in Fig. 1. From the figure, it is clear that the variations of 


these solutions are < 0.02% on binding energy and 0.01% on 
nuclear radii over the range of major shell Fermions Np from 
10 to 14. But this realative changes are reduced to < 0.002% 
on binding energy and 0.001% on nuclear radii for Np value 
from 14 to 22. Hence, the desired number of major shells for 
Fermions and bosons are Np = 18 and Np = 20 for the con¬ 
sidered mass region. However, the number of mesh points for 
Gauss-Hermite and Gauss-Lagurre integration are 20 and 24, 
respectively. For a given nucleus, the solution corresponding 
to maximum binding energy is treated as a ground state and 
other solutions are the intrinsic excited states of the nucleus. 
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TABLE III: Same as Table II, only for ^^®Mo and isotopes. 


N 

RMF (NL3) 

RMF (NL3*) 

Experiment 

RMF (NL3) 

RMF (NL3*) 

Experiment 


BE 

Teh 


BE rch 

/?2 

BE Vch 

P 2 

BE r^h 

P2 

BE 

Teh 

P2 

BE j32 

84-128MO 



88 

-132pd 


42 

696.7 

4.34 

-0.206 

696.2 4.34 

-0.203 



697.8 4.42 

0.002 

697.8 

4.42 

0.002 


44 

722.4 

4.34 

0.001 

722.6 4.34 

0.002 

725.8 


730.0 4.44 

0.094 

729.9 

4.44 

0.095 


46 

748.2 

4.33 

0.003 

748.2 4.33 

0.003 

750.1 


760.9 4.43 

0.101 

760.7 

4.43 

0.104 


48 

773.4 

4.33 

0.001 

773.1 4.33 

0.001 

773.7 


789.5 4.42 

0.004 

789.1 

4.42 

0.005 


50 

796.9 

4.33 

0.001 

796.4 4.33 

0.001 

796.5 4.3156 0.1058 

817.4 4.42 

0.001 

816.8 

4.41 

0.001 

815.0 

52 

812.7 

4.34 

0.001 

812.2 4.34 

0.001 

814.2 4.3518 0.1509 

836.9 4.43 

0.003 

836.3 

4.42 

0.004 

836.3 

54 

828.1 

4.38 

0.174 

827.1 4.37 

0.158 

830.8 4.3841 

0.1720 

855.9 4.46 

0.136 

855.4 

4.46 

0.139 

856.4 

56 

843.5 

4.42 

0.230 

842.2 4.41 

0.220 

846.2 4.4088 0.1683 

874.2 4.48 

0.176 

873.6 

4.48 

0.177 

875.3 4.4839 0.196 

58 

857.2 

4.45 

0.268 

855.7 4.43 

0.246 

860.5 4.4458 0.2309 

891.1 4.51 

0.189 

890.5 

4.50 

0.188 

892.8 4.5086 0.209 

60 

871.2 

4.50 

0.366 

869.1 4.49 

0.356 

873.9 

0.311 

906.9 4.52 

0.187 

906.2 

4.52 

0.184 

909.5 4.5322 0.229 

62 

883.6 

4.53 

0.386 

881.4 4.52 

0.382 

886.9 

0.362 

921.8 4.53 

0.190 

921.1 

4.53 

0.179 

925.2 4.5563 0.243 

64 

895.4 

4.49 

-0.234 

894.0 4.49 

-0.228 

898.9 

0.354 

936.1 4.57 

0.240 

934.9 

4.54 

0.165 

940.2 4.5776 0.257 


895.2 

4.54 

0.379 

893.9 4.54 

0.377 









66 

907.4 

4.51 

-0.236 

905.7 4.51 

-0.233 


0.38 

951.8 4.59 

-0.231 

950.4 

4.58 

-0.229 954.3 0.220 


906.5 

4.56 

0.374 

904.9 4.55 

0.373 



950.3 4.60 

0.292 

949.5 

4.60 

0.290 


68 

918.1 

4.53 

-0.241 

915.9 4.53 

-0.239 



965.7 4.62 

-0.23 

963.9 

4.60 

-0.234 967.6 0.164 









963.5 4.62 

0.304 

961.9 

4.60 

0.301 


70 

927.0 

4.55 

-0.231 

924.6 4.54 

-0.223 



977.9 4.63 

-0.221 

975.8 

4.61 

-0.226 

0.207 









975.3 4.60 

0.216 

974.1 

4.59 

0.216 


72 

935.7 

4.55 

-0.197 

933.2 4.55 

-0.190 



989.1 4.63 

-0.198 

986.8 

4.62 

-0.184 


74 

944.6 

4.57 

-0.180 

941.7 4.56 

-0.179 



1000.6 4.63 

-0.163 

998.1 

4.62 

-0.151 


76 

951.8 

4.58 

-0.172 

948.4 4.58 

-0.171 



1011.9 4.63 

0.115 

1009.8 

4.62 

0.114 


78 

958.4 

4.57 

0.114 

955.0 4.56 

0.094 



1022.8 4.64 

0.104 

1020.2 4.63 

0.105 


80 

965.1 

4.58 

0.042 

961.5 4.57 

0.037 



1032.4 4.65 

0.073 

1029.3 4.64 

0.076 


82 

972.2 

4.59 

0.002 

968.0 4.59 

0.001 



1041.5 4.65 

0.001 

1037.8 4.65 

0.002 


84 

973.7 

4.60 

0.002 

969.6 4.60 

0.029 



1044.5 4.66 

0.027 

1040.9 4.67 

0.043 


86 

975.8 

4.62 

0.114 

971.7 4.61 

0.101 



1048.7 4.69 

0.130 

1045.1 

4.69 

0.132 



A. Potential energy snrface 


Conventionally, in case of a quantum mechanical system, 
the path followed by the different solutions at various defor¬ 
mation define a potential barrier or potential energy surface, 
which can be used for the determination of the ground state 
of a nucleus. More elaborately, from the potential energy sur¬ 
face (PES) obtained from a self-consistent relativistic mean 
field theory, one can regulate the reasonable results for the 
ground state similar to the non-relativistic calculations GSl. 
Since quadrupole deformation plays the most important and 
dominant part, we have neglected the other deformation coor¬ 
dinates in the present study for simplicity and low computa¬ 
tion time cost. Here, the potential energy curve is calculated 
microscopically by the constrained RMF theory ifSTl l52l l79l - 


ISTl . The expectation value of the Hamiltonian llhTllTniSOll at 
certain deformation is given as, 


u' 


< > 


( 10 ) 


where A is the constraint multiplier and Hq is the Dirac mean 
field Hamiltonian. The convergence of the numerical solu¬ 
tions on the binding energy and the deformation are not very 
much sensitive to the deformation parameter /3o of the har¬ 
monic oscillator basis for the considered range due to the large 
basis. Thus the deformation parameter /?o of the harmonic os¬ 
cillator basis is chosen near the expected deformation to ob¬ 
tain high accuracy and less computation time period. 

The potential energy surface as a function of deformation 
parameter j32, for the proton rich nucleus ®^Zr, the double 
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FIG. 2: (Color online) The potential energy surface of * 2 , 90 , 110 , 1202 ^ 
as a function of quadrupole deformation parameter P 2 for both NL?> 
and NL3* forces in axially deformed relativistic mean field calcula¬ 
tions. 


FIG. 3: (Color online) The two neutron separation energy as a func¬ 
tion of neutron number from RMF theory with NL3 and NL3* force 
parameter for and nuclei 

are compared with the experimental data 1821 . 


magic nucleus and the neutron rich nucleus 110 T 202 J. ^re 
shown in Fig. 2, as a representative case. All other Mo—, 
Ru— and Pd— isotopes are also showing the similar behav¬ 
iors, which are not given here. The energy (Ei, = Eg ^ — E^.s 
on the F—axis is the difference between the ground state en¬ 
ergy to other constraint energy solutions. The solid and dotted 
line in the figure are for NL?> and NLi* force, respectively. 
The calculated PES for both the cases are shown for a wide 
range from oblate to prolate deformations. We notice from the 
figure that there are more than one minima appear at different 
P 2 ■ The magnitude of binding energy for the corresponding 
minima shows that the ground state solution appear at a cer¬ 
tain value of j32- The (32 for the ground state is not same for all 
isotopes of Zr (see Table. II-III. For example, the ground state 
solutions for ^^Zr, ^^Zr, iio.i20zr and i^Ozr are ~ -0.2, 0.0, 
0.4 and 0.0, respectively. One can find similar nature for both 
the force parameter, hence one can conclude that the ground 
state properties of these nuclei are independent of the force 
parameters used. 


B. Nuclear Binding energy and quadrupole deformatiou 

The calculations mainly explain the nuclear structure as 
well as the sub-structure properties, based on the basic in¬ 
gredients such as binding energy (Sb), quadrupole moment 
Q20, nucleonic density distribution p(rj^,z) = Pp{r 2 _,z) + 
Pn(f_L 5 z), and rms nuclear radii etc. Nevertheless, the present 
study demonstrates the applicability of RMF on the nuclear 
structure study for transition nuclei near neutron drip-line. 
The obtained results for binding energy BE, the quadrupole 
deformation parameter (32 and the charge radius Vch for NL3 
and NL3* force parameter for the isotopic chain of Zr, Mo, 
Ru and Pd are listed in Table II-III along with the experi¬ 
mental data ll82l . It is worth mentioning that the obtained 
results from the NL3 force parameter almost matches to the 


mass table by G. A. Lalazissis et al. ifSOl except few nuclei, 
though the small difference are acceptable in the accuracy of 
mean field level. Further, we notice on the binding energy and 
the rms Vch for all nuclei over the isotopic chain from RMF 
agree well with the experimental values. Quantitatively, the 
mean deviation of BE and Vch between the calculated result 
and the available experimental data over the isotopic chain 
are ^ 0.01 and 0.004, respectively. Further, the quadrupole 
deformation parameter (32, for both ground ig.s.) and selec¬ 
tive excited states (e.s.) are also given in Table II-III. In some 
of the earlier RMF and Skyrme Hartree-Fock (SHF) calcu¬ 
lations, it was shown that the quadrupole moment obtained 
from these theories reproduce the experimental data pretty 
well ll^l5ni52ll^l7ni74l[8TH85l. From the table, one can 
find that the shape of few nuclei is not consistent with the 
experimental observed shape. In this context, we have also 
estimated the first excited state solution for these nuclei cor¬ 
respond to the experimental deformations (see Table II-III). A 
careful inspection to these solutions shows that the small dif¬ 
ference in the binding energy is an indication of shape coexis¬ 
tence. In other words, the two solutions in these nuclei are al¬ 
most degenerate and might have large shape fluctuations. For 
example, in ®^Zr the two solutions for (32 = -0.197 and (32 = 
0.25 are completely degenerate with binding energies of 691.3 
and 691.0 MeV, respectively. Hence, the ground state can be 
changed to the excited state and vice verse by a small change 
in the input, like the pairing strength, etc., in the calculations. 
Similar behavior is also observed for few other nucleus in the 
present analysis are listed in Table II-III. Such phenomenon 
is known to exist in many other mass regions of the nuclear 
chart in [871 ■ 

C. Two neutron separation energy S2n (Z,N) 

Two neutron separation energy S 2 n (Z, N), can be estimated 
from the ground state nuclear binding energies of BE{Z, N), 
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FIG. 4: (Color online) The differential variation of the two neu¬ 
tron separation energy dS 2 n as a function of neutron number from 
RMF theory with NL3 and NL3* force parameter for 
84-128mo, 88-i32pj nuclei are compared with the 

experimental data l82l . 



FIG. 5: (Color online) The root-mean-square charge distribution rch 
of ®®"^“Ru and nuclei from RMF 

theory with NL3 and NL3* force parameter are compared with the 
experimental data 1891 . 


BE{Z, N — 2) and the neutron mass to„ with the relation: 

S 2 n{Z,N) = -BE{Z,N) + BE{Z,N-2) + 2mn, (11) 

The BE of the ^Xz and ^~^Xz are calculated from RMF 
for NL3 and NL3* force parameters. It is essential to have 
very precise mass measurements to predict the correct esti¬ 
mation of the nucleon separation energy S 2 n- The calculated 
S 2 n energy from RMF as a function of neutron number for 
Zr, Mo, Ru and Pd isotopes are compared with latest exper¬ 
imental data ll82l . shown in the Fig. 3. From the figure, it 
is clear that in an isotopic chain, the S 2 n energy shows the 
well-known regularities for a given atomic number i.e. the 
S 2 n decreases smoothly as the number of neutron increases 
in an isotopic chain. A sharp discontinuities (in other word 
kinks) appears at neutron magic numbers at = 50 and 82. 
In energy terminology, one can write, the energy necessary to 
remove two neutrons from a nucleus (Z, Nmagic+Z) is much 
smaller than that to remove two neutrons from the nucleus (Z, 
Nmagic), which breaks the regular trend. 


D. Differential variation of two neutron separation energy 


The differential variation of the two neutron separation 
energy (S' 2 n) with respect to the neutron number (N) i.e. 
dS 2 n{N, Z) is defined as 


dS2n{Z,N) 


S2n{Z,N E 2)- S2n{Z,N) 
2 


( 12 ) 


The dS 2 n{X,Z) is one of the key quantity to explore the 
rate of change of separation energy with respect to the neu¬ 
tron number in an isotopic chain. Here, we are calculated the 
dS 2 n{N, Z) for NL3 & NL3* force parameter. Further, we 
have also estimated the dS 2 n (N,Z) energy from the experi¬ 
mental S 2 n energy. In Fig. 4, we are compared the exper¬ 
imental values with our calculation for Zr, Mo, Ru and Pd 


isotopes. In general, the large sharp deep fall in the dS 2 n over 
an isotopic chain shows the signature of neutron shell closure. 
In other word, this deviation in the general trend may disclose 
some additional nuclear structure features. From the figure, 
we observed the same characteristics for all Z=38-46. 


E. The root-mean-square charge distributions 

The root mean square (rms) matter radius from relativistic 
mean field theory can be expressed as: 

{rD = jj p{r±,zydT, (13) 

where A is the mass number and p{r±,z) is the axially de¬ 
formed density. The rms charge radius can be calculated from 
the rms proton radius (r^) with simple algebraic relation, 

(rL) = {rl) + 0.64. (14) 

From the theoretical point of view, the macroscopic- 
microscopic models ll88l and microscopic mean-field formu¬ 
lations using effective interactions are most sophisticated ap¬ 
proaches to determine the rms charge radius in comparison 
with experimental data ll89ll . In this present work, we have 
shown the variations or fluctuations of the charge radii on the 
top of a fairly smooth average behavior in an isotopic chain. 
The results from RMF approaches for NL3 and NL3* pa¬ 
rameters along with the available experimental data are shown 
in Fig. 5. From the figure it is clear that the obtained radii 
from RMF for ®4-i28]v[o, 86-130 r^ 88-i32pj 

follows closely the experimental data lf89l . For most of the 
nuclei, the experimental values are unavailable, the RMF pre¬ 
diction are made for the charge radius of such a nucleus that 
awaits experimental confirmation. The circle, square and tri¬ 
angle symbols indicate the ground state states data for NL3, 
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FIG. 6: (Color online) The single particle energy levels ti for ®°Zr, 
^°'*Zr, ^^'‘Zr and ^'^^Zr from RMF model with NLS force parameter 


NL3* and experiment, respectively. Further, the solid circle 
and solid square symbols indicate the shapes correspond to 
the hrst intrinsic excited states obtained from NL3 and NL3* 
force. From the hgure, one can observe the smooth behavior 
for lighter isotopes, then there is a small fall in the charge radii 
for Zr, Mo, and Pd at about N=62, 64, 72 and 74. This fall 
corresponds to the transition from the prolate to the oblate and 
vice versa. But the magnitude for both the states are different, 
i.e., the oblate deformation is at /32 — 0.2 while the prolate 

one appears with f32 ~ 0.4. In case of Pb isotope, the change 
is the radii only at one place i.e. at A^=74. Further, one can no¬ 
tice that the tiny change in the calculation can lead to the first 
intrinsic excited state as a ground state (see Fig. 2). In other 
word, we can practically degenerate the ground state binding 
energy for the deformation corresponding to the first intrinsic 
excited state. Thus, the inconsistency in the could be ex¬ 
plained in terms of conhguration mixing i.e the actual ground 
state is not only the spherical configuration but also from the 
neighbor deformed intrinsic excited states. 


F. Single particle energy levels 

The above analysis, we found some signatures of shell clo¬ 
sures at N= 82 for all these isotopic chains. As a further con¬ 
firmatory test, the single-particle energy levels for neutrons 
in isotopic chains are examined. The obtained single particle 
levels ei for even isotopes of Zr near Fermi level are shown 
in Fig. 6 for NL3 force as an ideal case. However, we have 
obtained similar results for all isotopic chain from NL3 and 
NL3* force parameter. We observed the large gap at N = 82, 
near the drip-line region. In general, the spin-orbit splitting of 
the levels are scale down for neutron-rich nuclei, but lhii /2 
level (at N=82) is higher in the Zr nuclei studied. Quantita¬ 
tively, in 3222r, the Ae^ = ei(l*i3/2)ei(1^9/2) at N = 82 is 
4.5 MeV, which is a considerably large value compare to the 
neighbor splitting. Almost identical behavior is noticed for 
the isotopic chain of Mo, Ru and Pd nucleus at N=82, irre¬ 
spective of force parameter used. Such a rearrangement of the 


Neutron Proton 



X = Y 


FIG. 7: (Color online) The contour plot of the axially deformed 
ground state density distribution of proton and neutron for ®^Zr and 
«0Zr. 


single-particle orbitals at N=82, well accepted the shell clo¬ 
sure at N=82 for the considered transitional nuclei. 


G. The contour plot of the axially deformed densities 

In the above figures and tables, we have shown the results 
some of the structural observables such as binding energy, 
quadrupole deformation, root-mean-square radius, separation 
energy, differential separation energy and single-particle en¬ 
ergy levels in comparison with the experimental data 182] [89), 
wherever available. Here, we have focused on the ground 
of 84-128]y[Q^ 86-130]^^ 88-132p^ nuclei along 

with few selectively excited states. Based on these structural 
observables, we found some significant signature of the shell 
closure at A = 82 (drip-line region) in the isotopic chains. 
Further, the abnormal change in the S 2 n and dS 2 n in the iso¬ 
topic chain of Zr, Mo and Ru nucleus suggest a shape co¬ 
existence at A ~ 64 and 74. This divergence over an isotopic 
chain can be cut down by taking the dynamical correlations 
beyond mean-field ll90l - l^ . 

To get a complete picture into the reason behind such dis¬ 
crepancy over the isotopic chain, we have shown the contour 
plot of the axially deformed density of proton and neutron of 
these nuclei. In Fig. 7-8, we have displayed the distribution of 
Zr isotopes for N=42, 50, 60 and 82 as representative cases. 
All the isotopes of Mo, Ru, and Pd also showing similar be¬ 
havior as Zr as shown in Fig. 7-8. From the figure, one can 
clear identify the spherical, oblate, prolate shapes correspond¬ 
ing to their ^2 values as the local minima in the PECs. Similar 
calculations can also be found in Ref Il93]l94ll . In these fig- 












































9 


Neutron Proton 



X = Y 


FIG. 8: (Color online) The contour plot of the axially deformed 
ground state density distribution of proton and neutron for and 
i 26 Zr. 


ures, we can see that the transition from oblate to prolate at 
N=42, then change to the spherical structure at N= 50 and 
further changing the deformations to prolate one. Even the 
proton number is fixed in the isotopic chain, still we found 
a little change in the density distribution due to the influence 
of excess neutron number. Following the color code, the red 
and light gray color corresponding to the high density 0.09 
fm“^) and low density (~ 0.00 lfm“^), respectively. More in¬ 
spection on the figures shows that the central density of the 
proton increases as compared to the neutron with respect to 
the neutron number. In this region, few isotopes of Mo (for 
ii6 “1 i®Mo) are the triaxial shape in their ground state, which 


is very close to the axial solutions ll94ll . In other word, the 
location of minima for a triaxial solution for these isotopes of 
Mo are almost same as the minima appear for an axial prolate 
axial solution. Hence, we have used the simple axial deformed 
calculation, which is adequate for the a qualitative descrip¬ 
tions of structural observables in this mass region. 


IV. SUMMARY AND CONCLUSIONS 

We have used self-consistent relativistic mean held theory 
with most popular NL3 and recent NL3* force parameters to 
study the structural evolution in transition nuclei. The conjec¬ 
ture has been made from the binding energy, neutron separa¬ 
tion energies, differential variation of separation energy, the 
root-mean-square charge radii and the single particle energy 
levels of these nuclei. In this present calculations we have 
shown that Zr, Mo and Ru isotopes undergo a transition from 
oblate to prolate shapes at ~ 64 and 74. But, in case of Pd 
follows a smooth pattern through out the isotopic chain. We 
have also shown the dependence of nuclear charge radii on 
deformation also play an crucial role on their structural tran¬ 
sition. Further, we have also observed a large shell gap at = 
82 near drip-line region, almost same in magnitude at N= 50 
for these considered nuclei, which is a well-known feature for 
mean-field calculation. We have also demonstrated the effi¬ 
ciency of RMF theory calculations to reproduce those features 
and therefore to make predictions in unexplored regions. 
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